Improved di-neutron cluster model for He scattering 
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The structure of the three-body Borromean nucleus 6 He is approximated by a two-body di-neutron 
cluster model. The binding energy of the 2n-a system is determined to obtain a correct description 
of the 2n-a coordinate, as given by a realistic three-body model calculation. The model is applied to 
describe the break-up effects in elastic scattering of 6 He on several targets, for which experimental 
data exist. We show that an adequate description of the di-neutron-core degree of freedom permits 
a fairly accurate description of the elastic scattering of 6 He on different targets. 

PACS numbers: 24.10.-i, 24.10. Eq., 25.10.+S, 25.45. De, 25.60. Gc 



I. INTRODUCTION 

The scattering of a weakly bound projectile by a target 
represents a challenging as well as an interesting problem 
in nuclear physics. A proper understanding of the pro- 
cess requires an accurate description of the structure of 
the projectile, including all bound and unbound states 
that can be effectively coupled during the collision. In 
the case of weakly bound two-body projectiles the prob- 
lem has been solved using the Continuum Discretized 
Coupled-Channels (CDCC) method [3, 0,0]. Within the 
CDCC method, the reaction process of a loosely bound 
two-body projectile by a structureless target is treated 
within a three-body picture. The idea of the method is 
to represent the continuum part of the two-body projec- 
tile spectrum by a finite set of square integrable states. 
These states are then used to generate the diagonal as 
well as non-diagonal coupling potentials that enter the 
system of coupled equations. 

In principle, the method can be extended to three-body 
projectiles. This will be the case, for instance, of Bor- 
romean nuclei, consisting of three-body loosely bound 
and spatially extended systems, typically composed of a 
compact core plus two weakly bound neutrons (n + n + c), 
and with no bound binary subsystems. In this 
description of the three-body spectrum of the projectile 
is required. However, the calculation of the unbound 
spectrum of a three-body system is a very complicated 
problem by itself. In general, each physical state will 
be a complicated superposition of many channels with 
all possible spin and orbital angular momenta configura- 
tions. The calculation of the coupling potentials and the 
solution of the set of coupled equations in this large ba- 
sis represents a complicated task. Despite these difficul- 
ties, in two recent works @,0|, this method has been ap- 
plied to describe the scattering of 6 He on 12 C and 209 Bi. 
These calculations reproduce successfully existing elastic 
scattering data for these reactions and represent an im- 
portant advance towards the understanding of few-body 
nuclear reactions. 

Most of the complexity of these processes involving 



Borromean nuclei arises from the fact that these sys- 
tems exhibit many excitation modes, which can be as- 
sociated with two different degrees of freedom: the n — n 
relative motion, and the (nn) — c motion. In general, 
both modes will be excited during the collision. How- 
ever, when the system is scattered by a medium mass 
or heavy target, the projectile-target interaction will ex- 
cite mainly the coordinate between the neutrons and the 
core, since the repulsive Coulomb interaction will tend to 
repel the charged core, while the neutrons can approach 
closer to the target. Moreover, the nuclear interaction 
will attract more strongly the weakly bound neutrons. 
So, the net effect of the interaction with the target will 
be to stretch the nn—c coordinate, pushing the core apart 
from the target and pulling the neutrons close to it. Thus, 
a description of the projectile excitation mechanism that 
takes into account explicitly the nn — c coordinate should 
explain the main features of the reaction mechanism of 
the three-body system with the target. 

Given the complexity of the full CDCC calculations 
with three-body projectiles, the development of these 
simple models can be very helpful to understand the main 
features of these processes by retaining only the essential 
ingredients to keep the model realistic. 

In this work we revisit the so called di-neutron model 
for the 6 He case. In Sec. [HI we address the problem of 
the 6 He structure within a three-body model. In Sec. lIII[ 
we review the di-neutron model for this nucleus, and we 
propose a method to improve the accuracy of the model, 
while keeping its simplicity. In Sec. |IV] the new method 
is tested against existing experimental data for 6 He scat- 
tering on several targets. Finally, Sec. V is for summary 
and conclusions. 

II. THREE-BODY MODEL FOR 6 HE 

Within a three-body picture, the wavefunction of the 
6 He system can be conveniently expressed in terms of 
one of the Jacobi sets of coordinates. For the purposes 
of the present work, the most suitable representation is 
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that in terms of the neutron-neutron relative coordinate, 
x, and the nn- 4 He coordinate, y. This wavefunction, 
here denoted 4 ,3B (x, y), can be obtained by solving the 
Schrodinger equation, using any of the methods proposed 
in the literature. Here, we followed the procedure pro- 
posed in [f| @|, in which the wavefunction is expanded 
in hyperspherical coordinates. The basic ingredient of 
the calculation are the two-body interactions between 
the subsystems (n — n and n — a). Besides the two- 
body potentials, the model Hamiltonian also includes a 
simple central three-body force depending on the hyper- 
radius. This is introduced to overcome the under-binding 
caused by the other closed channels, such as the t+t 
channel. The n- 4 He potential is taken from Refs. [1, @|, 
with central and spin-orbit components, and the neutron- 
neutron potential, with central, spin-orbit and tensor 
components, from the prescription of Gogny, Pires and 
Tourreil [XQ] . These calculations were performed with 
the code STURMXX [Tl| which uses the formalism de- 
scribed in Q. The maximum hyperangular momentum 
was set to i^ max = 20 and the three-body force was 
adjusted to give the right binding energy. The calcu- 
lated three-body wave function has a binding energy of 
£6 = 0.955 MeV and a rms point nucleon matter radius 
of 2.557 fm when assuming an alpha-particle rms matter 
radius of 1.47 fm. Further details of these calculations 
can be found in Ref . 0, [12] • It should be noted that the 
three-body wavefunction is a complicated superposition 
of many channel configurations. Each channel is charac- 
terized by the angular momentum in the n-n and (nn)-a 
coordinates (l x and l y ), the total orbital angular momen- 
tum (L) and the total spin of the neutron pair (S x ). In 
the 6 He ground state the dominant configuration corre- 
sponds to S x = l x — l y =0, which contributes to 80% of 
the norm. 

In order to compare with the di-neutron model, pre- 
sented below, we consider now the behavior of the wave- 
function in the y coordinate. For this purpose, we calcu- 
late the probability density in the (nn)-a relative coor- 
dinate, y, here denoted as p(y). This was calculated by 
integrating the square of the three-body wavefunction on 
the neutron-neutron coordinate, i.e. 



P ZB {y) = y z / |* aa (x, y )|^xdfi 



3B, 



(1) 



where ^(x, y) is the total three-body wavefunction and 
fly denotes the angular variables (0 y , 4> y ). This density is 
plotted in Fig.Q]with the thick solid line. To illustrate the 
dominance of the S x = lx = ly = component, we show 
also the same quantity, retaining only this component in 
the wavefunction (thick dashed line). It can be noticed 
that, for y > 5 fm, the di-neutron density is completely 
determined by this component. Consequently, a realis- 
tic model for the 6 He ground state wavefunction must 
account, at least, for this configuration. 

In scattering calculations involving the 6 He nucleus it 
is essential to include also a realistic description of the 
continuum states, given the large breakup probability of 
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PIG. 1: (Color online) Neutron density in the di-neutron 
model, compared with a realistic three-body model. The di- 
neutron calculations use a Woods-Saxon potential with radius 
R = 1.9 fm and diffuseness a = 0.25 fm. 



weakly bound nuclei. In the case of the Coulomb inter- 
action, the response of the continuum to excitations of 
multipolarity A is conveniently treated in terms of the 
reduced transition probability, B(EX) [13]. In Fig. [2] we 
consider the B(E1) (upper panel) and B(E2) (bottom 
panel) distributions, plotted as a function of the excita- 
tion energy of the 6 He nucleus with respect to the ground 
state. In both panels the full three-body calculation is 
depicted by the thick solid line. In these calculations, 
the continuum states were represented by true scattering 
wave functions, as reported in Ref. [§]. The narrow peak 
in the B(E2) corresponds to the known 2 + low lying res- 
onance. 



III. THE DI-NEUTRON CLUSTER MODEL OF 
6 HE 

We want to consider situations in which the (nn) — c 
degree of freedom is more relevant than the nn degree of 
freedom in 6 He. This is the case, for example, when elec- 
tric operators are considered in structure calculations, or 
when Coulomb forces dominate in a collision. One could 
think then in approximating the three-body wavefunc- 
tion by a product of two-body wavefunctions, i.e. 



* 3B (x,y)^^(y)^(x) 



2B, 



(2) 



The di-neutron model takes into account only the 
(nn) — c degree of freedom, whereas the relative mo- 
tion between the two valence neutrons is ignored. This 
amounts to consider that the neutron pair remains in 
a highly correlated state "0( x ) during the collision and, 
hence, excitations in this coordinate are not permitted. 
Moreover, the neutrons are assumed to be coupled to spin 
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FIG. 2: (Color online) B(E1) and B(.E2) transition strengths 
for the 6 He nucleus using different structure models. The solid 
lines correspond to the three-body calculation obtained with 
the true scattering states [§]. Di- neutron calculations using 
for the 2n-a binding energy either the two-neutron separation 
energy (dotted-line) or the modified value et = —1.6 MeV 
(dotted-dashed line) are also shown. 



zero, and bound to an inert a core in a s-wave relative 
motion which, as we have seen before, is the dominant 
configuration in the 6 He ground state wavefunction. This 
model is inspired in the deuteron-a cluster model, which 
has been very successful in the description of 6 Li scatter- 
ing within the continuum discretized coupled channels 
method (eg. 0, [H, 13). 

The problem in the di-neutron model lies on evalu- 
ating the wavefunctions, for the bound and continuum 
states, describing the motion of the halo neutrons rela- 
tive to the a core, ^f 2B (y). Following the cluster model, 
one assumes that these wavefunctions can be obtained 
as the eigenstates of the Hamiltonian corresponding to 
a certain 2n-a interaction. Typically, one assumes some 
reasonable geometry for the 2n-a interaction, and then 
adjusts the potential depth to obtain a given binding en- 
ergy for the 2n-a system. 

In all the calculations here presented, the a + 2n in- 
teraction was parameterized using a standard Woods- 
Saxon form, with radius i?o=l-90 fm and diffuseness 



a = 0.25 fm, which corresponds to the set III of Ref. [171 ] . 
The ground state wavefunction was assumed to be a pure 
2S configuration, since, due to the Pauli principle, the IS 
state is forbidden. 

So, the key question is: which is the binding energy 
that one should use for the 2n-a system, so that the cor- 
responding wavefunction gives a reasonable description 
of 6 He in a di-neutron model? 

In the a ppli cation of the deuteron-a cluster model to 
6 Li (l3.[lH.Tl3|. one evaluates the binding energy for d — a 
just as the separation energy of 6 Li into d + a. This is 
a reasonable procedure, which can be applied because 
the deuteron is bound by 2.2 MeV, which is more than 
the separation energy of 6 Li into d + a, and so one can 
argue that the relative wavefunction of the valence proton 
and neutron within 6 Li is not very different from that in 
a free deuteron. Besides, the deuteron-a cluster model 
gives reasonable values for the mean square radii of 6 Li. 

In the a ppl ications of the di-neutron model done so 
far to 6 He [H [TJ, O, S3, El HI, the binding energy of 
the di-neutron has been taken as the two-neutron sepa- 
ration energy of 6 He, i.e., |ef,| =52^=0. 975 MeV. It should 
be noticed that, with this choice, one is assuming implic- 
itly that the relative wavefunction of two neutrons within 
6 He would be in a state similar to that of two neutrons 
with zero relative energy. This leads to an unrealistic 
wavefunction for the di-neutron-cv motion, as discussed 
below. 

The use of a binding energy ej, = —0.975 MeV yields 
the potential depth V = 93.51 MeV. For the i = and 
1=1 continuum states we used the same potential as for 
the ground state. For the i = 2 continuum, the potential 
depth was changed to Vb = 91.25 MeV, in order to get the 
2+ resonance at the value obtained in the the three-body 
calculation which, in turn, is close to the experimental 
value (e x = 0.825 MeV above the breakup threshold). 

To illustrate this, we compare the density probability 
associated with the nn — c coordinate in the two- and 
three-body models. In the di-neutron model, the neu- 
tron density analogous to Eq. {1} is simply obtained as 
P 2B {y) — |y^(j/)| 2 ) where R{y) is the radial part of the 
wave function ^> 2B . In Fig. [TJ the density probability 
obtained with this model is given by the dotted line. 
When compared to the realistic three-body calculation 
(thick solid line) it becomes apparent that the former ex- 
tends to considerably larger distances. For example, the 
rms associated to the di-neutron-a coordinate is 4.36 fm, 
considerably larger than the prediction of the three-body 
model, 3.25 fm. In view of this result it is not surprising 
that the coupling of the ground state wavefunction with 
the dipole continuum states is unphysically enhanced in 
the two-body model. 

This is indeed the case, as we can see in Fig. [21 In 
both panels, the dotted line corresponds to the di-neutron 
model. These distributions clearly overestimate both the 
El and E2 strengths predicted by the three-body model 
(thick solid lines). Not surprisingly, previous attempts 
to describe 6 He scattering by heavy targets [23l . [2~3 | us- 
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ing this model showed that this simplified description of 
the 6 He nucleus tends to overestimate the effect of the 
continuum couplings. In [22] it was shown that, reducing 
the strength of the dipole couplings, the agreement with 
the data could be significantly improved. 



We consider that the 2n-a binding energy used in the 
di-neutron model should not be given by the two-neutron 
separation energy. The di-neutron system which appears 
in 6 He is a correlated state, which, in the absence of the 
a particle, will be given by a wave packet with positive 
expectation value of the energy. Thus, the actual 2n-a 
binding energy should be more negative, to compensate 
for the positive energy of the di-neutron. We propose 
to obtain 2n-a binding energy to reproduce the known 
properties of the 6 He system, such as the rms radius and 
the transition strengths, within the di-neutron model. 
Asymptotically, the di- neutron wavefunction behaves as 
cx exp(— ky) with k = y/2fi\eb\/h and fi the 2n-a reduced 
mass. Then, a natural choice for |e&| is to make the slope 
as close as possible to the three-body case. This leads 
to the value |e&| = 1.6 MeV. The density calculated with 
this value, shown by the dotted-dashed line in Fig. [TJ 
reproduces very well the three-body calculation (thick 
solid line) for separations beyond 4 fm. The di-neutron-a 
mean square separation obtained with the new wavefunc- 
tion is reduced to 3.4 fm, in much better agreement with 
the three-body result. With the new binding energy, the 
depth of the s-wave is modified to Vb=96.06 MeV. Again, 
this depth was used for the p-waves. The depth of the 
1 = 2 potential had to be changed to Vo=92.7 MeV, in 
order to get the 2 + resonance at the correct excitation 
energy with respect to the ground state. 



We note that, among the three geometries proposed in 
Ref. namely, the set I (a = 0.65 fm), set II (a = 
0.39 fm) and III (a = 0.25 fm) the latter is found to 
reproduce more accurately the three-body density. The 
other two geometries, having a larger diffuseness, give 
rise to a higher rms, even after modification of the two- 
neutron separation energy to correct the slope of the di- 
neutron density. 



It can be seen that, with this geometry and binding 
energy, the El and E2 transition strengths are also well 
reproduced. This is shown in Fig.[2]by the dotted-dashed 
lines. It is observed that this increase of the binding en- 
ergy reduces both strengths showing a much better agree- 
ment with the prediction of the three-body calculation. 
Note that, due to the modification of the binding energy, 
the breakup threshold appears at a higher energy in our 
di-neutron model. So, the di-neutron model does not de- 
scribe the low-energy continuum of 6 He, which is below 
1.6 MeV excitation energy. However, it describes fairly 
well the continuum around the maximum of the B(E1) 
distribution (2 MeV) and beyond. 
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FIG. 3: (Color online) Elastic scattering angular distribu- 
tion, divided by the Rutherford cross section, for the re- 



action °He+ 64 Zn at E 



10 and 13.6 MeV. The dashed 



line is the cluster-folding calculation without inclusion of the 
continuum. The dashed line and the thick solid line are 
the di-neutron model calculations using the binding energy 
£b = Sm = —0.975 MeV, and the modified binding energy 
£6 = — 1.6 MeV, respectively. The experimental data are taken 
from El. 



IV. CALCULATIONS 

In the remaining, we compare the two-body models 
discussed above with the elastic scattering data for sev- 
eral reactions induced by 6 He. In all the calculations 
here presented, we used the geometry of the 2n-a poten- 
tial with the smaller diffuseness (a = 0.25 fm). All these 
calculations are performed within the standard CDCC 
method Q. 

We first consider the reaction 6 He+ 64 Zn at Coulomb 
barrier energies, which was recently measured by Di 
Pietro et al. [25|]. The 6 He (=2n + a) continuum was 
discretized into N — 7 energy bins, evenly spaced in the 
asymptotic momentum k, and up to a maximum excita- 
tion energy of e max = 7 MeV. We included s, p and d 
waves for the 2n-a relative orbital angular momentum. 
Inclusion of / waves had a negligible effect on the elastic 



angular distributions. 

In these calculations, the a+ 64 Zn interaction was 
taken from the optical model fit performed in j25|. For 
the 2n+ 64 Zn interaction, we used the parameters of the 
d+ 56 Fe potential obtained in Ref. [26]. Diagonal as 
well as non-diagonal potentials were derived from these 
potentials by means of a single-folding method, as de- 
scribed elsewhere The coupled equations were in- 
tegrated up to 100 fm, and using 50 partial waves for 
the projectile-target relative motion. These calculations 
were performed with the code FRESCO |27l |. 

In Fig. [3] we show the results for the elastic scattering 
angular distribution at the laboratory energies E=10 and 
13.6 MeV, along with the data of Di Pietro ei al. For 
each energy three curves are shown: the dashed line is the 
cluster-folded calculation in which the projectile-target 
interaction is folded with the ground-state density of the 
6 He nucleus, without inclusion of the continuum. At 
the higher energy this calculation exhibits a pronounced 
rainbow which is not observed in the data. Moreover, 
at both energies these calculations clearly underestimate 
the data at backward angles. Inclusion of the continuum 
within the conventional di-neutron model (dotted lines) 
improves the agreement at backward angles, but reduces 
too much the cross section at the rainbow, thus under- 
predicting the data. This effect is a direct consequence 
of the overprediction of the B(E1) distribution in the 
di-neutron model, as explained above. Finally, the thick 
solid line is the CDCC calculation with the di-neutron 
model with a modified binding energy (eb = — 1.6 MeV). 
This calculation improves the agreement at the rainbow, 
particularly at E = 13.6 MeV. At E = 10 MeV this 
calculation slightly underestimates the data at backward 
angles, but we could not find an explanation for this dis- 
crepancy. 

Next, we study the 6 He+ 208 Pb^ 09 Bi, 197 Au reactions, 
which were recently analyzed in [22]. In Fig. [4] we present 
the calculations for the Pb and Bi targets, along with the 
experimental data from Kakuee et al. [28| and Aguilera 
et al. [29], respectively. Details of the fragment-target 
optical potentials and binning scheme can be found in 
Ref. [22|. The meaning of the lines is the same as in 
Fig.[3l In both cases, the calculation without continuum 
displays a marked rainbow, which is not observed in the 
data. We have added also a calculation without contin- 
uum, but with the modified binding energy et = —1.6 
MeV (thin solid line). This calculation illustrates the 
static effect caused by the change in the ground state 
wavefunction produced by the modification of the bind- 
ing energy. Qualitatively, these two calculations are very 
similar. In particular, the pronounced rainbow is still 
present in the new calculation. The similitude between 
these two calculations indicates that the trend of the 
data can not be simply explained by changing the size 
of the di-neutron— a wavefunction and, consequently, dy- 
namical effects are indeed very important. Inclusion of 
the continuum within the conventional di-neutron model 
(dotted line), produces a strong reduction of the cross 
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part). The dashed line is the calculation without inclusion of 
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sertion at intermediate angles, largely underestimating 
the data. As noted above, this is caused by the overesti- 
mation of the dipole couplings in this model. In the mod- 
ified di-neutron model the rainbow is also suppressed, but 
the final result is in very good agreement with the data. 

Finally, we discuss the results for the 6 He+ 197 Au re- 
action at E = 27, 29 and 40 MeV, and compare with 
the data of [30]. The lines have the same meaning as in 
figures [3] and [H Similarly to the case of the lead tar- 
get, at E — 27, 29 MeV the one channel calculation ex- 
hibits a pronounced rainbow, which is almost absent in 
the data. This effect is very well accounted for in the 
modified di-neutron calculation. At E — 40 MeV, the 
rainbow is suppressed somewhat, but not completely, in 
the full CDCC calculation. The lack of data at the rele- 
vant angles does not permit to make strong conclusions 
about the existence of the rainbow, but the agreement be- 
tween the data and the calculation is fairly good where 
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FIG. 5: Elastic scattering angular distribution, divided by 
the Rutherford cross section, for the reaction 6 He+ 197 Au at 
E = 27, 29 and 40 MeV. The dashed line is the cluster-folding 
calculation without inclusion of the continuum. The dotted 
line and the thick solid line are the calculations using the di- 
neutron model with et = — S2n = —0.975 MeV and £(, = —1.6 
MeV, respectively. The experimental data are taken from [3(il |. 



the comparison is possible. 

It is interesting to note that the underestimation of the 
data in the conventional di-neutron model is more pro- 
nounced at lower energies. This is because at lower en- 
ergies dipole Coulomb couplings become more important 
and, as we showed before, these couplings are unphysi- 
cally enhanced in the conventional di-neutron model. 

All these calculations show that the proposed model 
describes fairly well the elastic data for different targets 
and could even be used as a predictive tool for reactions 
for which data do not exist. The good agreement with 
the data clearly supports the idea, anticipated in the in- 
troduction, that the degree of freedom which enters ac- 
tively in the elastic scattering of Borromean systems on 
medium mass and heavy targets, is that for the relative 



motion between the halo neutrons and the core. 
V. SUMMARY AND CONCLUSIONS 

In this work, we have studied the application of the 
di-neutron model to describe 6 He structure and scatter- 
ing. We find that, when the di-neutron model is applied 
assuming for the 2n-a binding energy the two-neutron 
separation energy of 6 He, the description of the structure 
of 6 He obtained is not in agreement with the results of a 
realistic three-body calculation. One obtains an unreal- 
istically long tail of the 2n-a relative wavefunction, and 
too large values of the B(E1) and B(E2) distributions. 
When this model of 6 He is used in scattering calcula- 
tions, the couplings between the ground state and the 
continuum states are overestimated, and this produces 
too much absorption from the elastic channel. 

We have proposed a modified di-neutron model, in 
which the 2n-a binding energy is set to reproduce the 
density distribution in the 2n-a coordinate given in a re- 
alistic 3-body model. We find that a 2n-a binding energy 
of 1.6 MeV produces a rms and B(El) and B{E2) distri- 
butions which are similar to those obtained in a realistic 
three-body calculation. 

The model has been tested for several reactions in- 
duced by 6 He, providing in all cases a very good descrip- 
tion of the elastic scattering data. These results indicate 
that, despite its simplicity, the model can provide a use- 
ful and reliable description of reactions involving the 6 He 
nucleus. Using an identical procedure, the method could 
be also extended to other Borromean systems, such as 
n Li or 14 Be. 

We would like to emphasize that the present model 
is not intended to replace the realistic three-body cal- 
culations for the scattering of Borromean nuclei. The 
development of these models, although numerically more 
demanding, are of great importance for a full quantita- 
tive understanding of these processes. However, we be- 
lieve that simple models, as those discussed here, are also 
very useful to provide us with a transparent physical in- 
terpretation of these collisions. 
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